Search for the flavour-changing neutral current decay 
D° — > in pp collisions at y/s = 7 TeV with CMS 



CM 



Daniele Pedrini 



ON BEHALF OF THE CMS COLLABORATION 

INFN Sezione di Milano-Bicocca, 1-20126 Milano, ITALY 

< 

V 

^ A search for the flavour-changing neutral current decay D° — > /jL + /j,~ is 

performed in pp collisions at y/s = 7 TeV using 90 pb _1 of data collected 

by the CMS experiment at the LHC. No evidence is found for this decay 

mode. The upper limit on the branching fraction B(D° — > jj + yT) is 

00 5.4 x 10" 7 at the 90% confidence level. 
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1 Introduction 



One promising way to search for physics beyond the Standard Model (SM) is to search 
for decay modes that are extremely rare or forbidden. The observation of these modes 
at rates exceeding the prediction of the SM could open a window onto New Physics 
(NP). The Flavour-Changing Neutral Current (FCNC) decays are rare decays which 
proceed via an internal quark loop in the SM but are forbidden at the tree level. In 
the SM, the FCNC decay D° — > n + ji~ is highly suppressed by the Glashow-Iliopolus- 
Maiani (GIM) mechanism and by a factor of m^/m 2 -^ due to helicity. This decay 
proceed via a VKbox diagram, which also contributes to D° — D° mixing. Theoretical 
estimates of this branching ratio are approximately 10~ 18 from short range processes, 
increasing to 10~ 13 when long-distance processes are included pQ. However NP models 
can enhance these estimates by several orders of magnitude [2J- This is why these 
decays are so attractive: any detection, given current sensitivities, will be a clear sign 
of NP. Large Hadron Collider (LHC) experiments have the possibility to detect these 
rare decay modes. Furthermore, as charm is an up-type quark, the search for FCNC 
in the charm sector is complementary to B and K decay searches. 

In this note, a search for D° — > fi + n~ is presented using a data sample of pp 
collisions at y/s = 7 TeV, corresponding to an integrated luminosity of ~ 90 _1 , 
collected by the Compact Muon Solenoid (CMS) experiment at the LHC. 



2 Analysis 

The strategy of this analysis is to measure the ratio of branching fractions, ^+^^0(^-^+1^+ ; 
in such a way that most of the systematic uncertainties cancel out (throughout this 
note charge conjugate state is implied). The main challenge of this analysis is to 
reconstruct the normalization mode, which has a much smaller trigger efficiency A 
feature of heavy flavour events is to have a secondary vertex separated from the 
primary vertex. For this reason the semileptonic and the dimuon analyses are essen- 
tially topological: they are based on a search for primary and secondary vertices in 
the event. In the case when multiple primary vertices are reconstructed, the vertex 
with the highest sum of p\ of the tracks in the vertex is selected, where pr is the trans- 
verse momentum. This is the default CMS choice. Events without a primary vertex 
or events in which the selected primary vertex has a x 2 probability less than 1% are 
discarded. A detailed description of the CMS detector can be found in Ref. [3]. The 
main subdetectors used in this analysis to reconstruct the topological configuration 
of the event are the silicon tracker (composed of pixel and strip layers) and the muon 
stations. The trigger plays an essential role. The events selected for this analysis are 
those with a muon having a px greater than a certain threshold, which varies with 
running conditions. The events are selected with a two-level trigger system. The 
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first level requires a good quality muon, while the high level trigger (HLT), using 
additional information from the silicon tracker, imposes a cut on the px of the muon 
(for example HLT_Mu3 means an event with a muon having > 3 GeV/ c). As the 
luminosity of the LHC increased the triggers were prescaled. During the 2010 data 
taking, the lowest pt unprescaled trigger was varied six times from 3 (HLT_Mu3) to 
15 GeV/c (HLT_Mul5). In the 2011, the HLT_Mul5 trigger remained unprescaled 
for the first 54 pb _1 of collected luminosity. Later CMS trigger configurations are 
too inefficient to reconstruct the normalization decay mode D° — > K~ pL + v. The data 
sample considered is therefore divided into seven run periods: six for the 2010 data 
taking and one for 2011. 

The D° — > ji + fjT analysis begins with two opposite sign muon candidates. Muon 
candidates are required to be reconstructed both in the silicon tracker and in the 
muon stations. The global fit must have a x 2 P er degree of freedom (DOF) less than 
3, to have a distance of closest approach to the primary vertex in the transverse 
plane less than 2 mm, and to be within the pseudorapidity region \r]\ < 2.1, where 
i] = — In [tan(6'/2)] and 9 is the polar angle with respect to the counterclockwise beam 
direction. In addition, one of the two muon candidates must match the muon which 
fires the trigger {trigger muon); this implies the following cuts for the different 

2010 run periods: px > 3, 5, 7, 9, 11 and 15 GeV/c. The 7 th run period is the 

2011 data taking with px > 15 GeV/c. The other muon (second muon) must satisfy 
Pt > 3 GeV/ c. The two muon candidates must form a vertex with at least a 1% x 2 
probability (CL se c > 1%)- If a good secondary vertex is found, the position of the 
primary vertex is recomputed, excluding the two muons from the vertex. A point- 
back cut is applied: the angle between the D° momentum vector and the line of 
flight of the D° (direction between the primary and secondary vertices) must satisfy 
cos a > 0.99. Finally the D° candidate is combined with a track (px > 0.6 GeV/c), 
which is given the pion mass and must originate from the primary vertex, to form a 
D* + candidate. Combinations with AM = M(fi + fi~n + ) — M(fi + fi~) exceeding 180 
MeV are discarded; if more than one candidate is found with AM < 180 MeV/c 2 , 
the candidate with AM closest to the nominal PDG jl] mass difference is chosen. 

The semileptonic decay mode reconstruction, developed by E691 [5J, is based 
on the decay chain D* + — > D°(K~ ^ + u)7i + whereby the v momentum can be re- 
constructed provided the D° direction (vector between the primary and secondary 
vertices) is sufficiently precise. With this technique the longitudinal component of 
the neutrino momentum (p z v ) can be determined and the energy and momentum of 
the (Kfiu) candidate can be calculated. This (Kfiu) candidate is then combined with 
one reconstructed track from the primary vertex to determine the D* + mass. The sign 
of this candidate pion must be equal to that of the muon. To arbitrate between the 
two possible (pi) solutions (two-fold ambiguity) one chooses the solution which gives 
the smallest D* + mass; Monte Carlo (MC) studies show this to be correct ~80% of 
the time. Finally, the mass difference, AM = M(D* + ) - M(D°), is determined. The 
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D° — > K~fi + v analysis begins by considering one kaon candidate (that is a track as- 
sumed to be a kaon) and an opposite sign muon candidate. The muon candidate must 
satisfy the same cuts used for the trigger muon of the D° — > fi + fi~ analysis, while 
the kaon candidate is a track reconstructed in the silicon tracker with x 2 /DOF < 3, 
\t]\ < 2.1 and pr > 0.8 GeV/c. As in the D° — > fi + fi~ analysis, the kaon and the 
muon candidates are required to form a secondary vertex with CL S ec > 1% and the 
position of the primary vertex of the event is recomputed excluding the kaon and the 
muon if they belong to the primary. Once the direction between the primary and the 
secondary vertex is known, the E691 technique described above is used to determine 
the v momentum. To form a D* + , if more than one pion candidate (that is a track 
originating from the primary with px > 0.6 GeV/ c) is found within the AM range 
(AM < 180 MeV/c 2 ), that with AM closest to the nominal PDG mass difference 
is chosen. Due to the charge relation between the /i and the ir of the decay chain 
D* + — > D°(K~ fi + u)7i + , each D° candidate selects a Right Sign (RS) D* candidate 
with K~fi + 7i + and a Wrong Sign (WS) D* candidate with K~fi + 7i~ . 

To reduce prompt background, the separation between primary and secondary 
vertices was required to be greater than three times the uncertainty on the separation 
{L/ai > 3). This cut value was optimized on the normalization mode D° — > K~[i + i> 
to reduce bias. 

In the WS sample there is no evidence of a signal, as expected. This WS sam- 
ple is used to model the background of the RS sample. Figure [T] shows AM = 
M(K-fi + un + ) - M(K-fi + u) for the D* + ->• D°(K-fi + u)n + analysis in the 7 th pe- 
riod (the period with the largest dataset). The superimposed fit to the unbinned RS 
data includes two Gaussian functions (the wider one is bifurcated to take into account 
the threshold on AM) with the same mean plus a background function obtained from 
the WS sample. The fit returns 16 458 ± 204 D° — > K~ /i + z/ candidates. 




D n 1 1 1 1 i 

140 145 150 155 160 165 170 175 
AM - M(Knvn) - M(Ku.v) (MeV) 



Figure 1: AM = M(K~ fi + u7i + ) — M(K~ for the RS sample with fit superim- 
posed. The points are the data, the blue line shows the fit function composed of 
two Gaussians plus the background function modelled by the WS. This background 
function is shown as a dashed red line. 
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For the D° — » analysis , with the additional cut |AM — AMpoal < 3 MeV 

where 3 MeV corresponds to ~ 3.6 times the mass resolution measured in MC, we 
obtain the invariant mass shown in Figure [2] (data from all seven periods). In 

more detail, Figure [3] shows the invariant mass for each trigger period. There 

is no evidence of the D° — > decay. 



CMS,\/s = 7TeV 




Dimuon invariant mass (GeV) 

Figure 2: The fi + n~ invariant mass distribution. The light red band shows the signal 
region, while the two light blue bands show the sidebands. 



The upper limit on the branching fraction D° — > fi + fi is determined using the 
following formula: 

B(D° -+ » V) < B(D° -+ K-^u) x x x x *rec(*>) 

N(Knu) a{fj,fi) e trig (/i/i) e T ec{M) 

(1) 

where B(D° ->■ K~n+u) = (3.30 ± 0.13) x 10" 2 (PDG) 0] is the normalization 
branching fraction, N(fifi) is the 90% CL upper limit on the D° — > fi + fx~ yield, 
N(KfMv) is the number of D° — > K~ ji + u candidates, and a and e are the acceptance 
and the efficiencies of the two modes. 

The Monte Carlo simulation is used to determine the acceptance and efficiency 
ratios for the signal and normalization mode. The MC event samples (one for each 
period) are generated with Pythia 6.409 [S], the unstable particles are decayed with 
EvtGen [7], and the detector response is simulated with Geant4 [S]. 

The seven periods are simulated with the corresponding triggers and run condi- 
tions, including simulation of the pile-up. The number of MC events in each period 
is proportional to the corresponding statistics of the data. The average number of 
reconstructed primary vertices ranges from 1.7 to 5.5 and is well matched by the 
simulation. 

The acceptance is defined by a = N acc /Ng e n, where iVacc includes events in 
which both tracks (K~fi + or /i~yU + ) have \r]\ < 2.1 at the generation level, while 
the denominator iVgen is the total number of signal decays generated. The ratio of 
acceptance for the two modes is a(K ' \w) / 'a(fifi) = 0.995 ± 0.006 (stat.). 
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Figure 3: The invariant mass distribution for each trigger period. The light 

red band shows the signal region, while the two light blue bands show the sidebands. 
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The trigger efficiency is defined by e^ r jg = N^g/N acc , where A^ r jg is the num- 
ber of events that pass a particular trigger. Table shows the ratio of the trigger 
efficiencies for D° — > K~n + v and D° — > fi + fi~ for each of the seven periods. As 
charm is produced at low p?, the trigger efficiencies are very low, especially for the 
normalization mode. 

The reconstruction efficiency is defined by erec = ^rec/^trig' wnere ^rec is the 
reconstructed signal yield after all the analysis cuts. Table [T] shows the ratio of the 
reconstruction efficiencies for D° — > K~ /i + z/ and D° — > n + fi~ in the 7 periods con- 
sidered. The ratio erec(-^A iZ/ )/ e rec(A t A t ) ) ranging from 2.2 to 0.8, is highly dependent 
on the trigger. 

Table 1: Ratio of trigger and reconstruction efficiencies in the seven periods used in 
the analysis, uncertainties are statistical only. 



Trigger 


etri g (^^)/etrig(/ i / i ) 


e Tec (K nv) 1 e r ec(^) 


HLT_Mu3 

HLT_Mu5 

HLT_Mu7 

HLT_Mu9 

HLT_Mull 

HLT_Mul5 (2010) 

HLT_Mul5 (2011) 


0.149 ±0.002 
0.112 ±0.002 
0.102 ±0.003 
0.099 ±0.003 
0.097 ±0.003 
0.085 ±0.003 
0.087 ±0.002 


2.215 ±0.197 
1.651 ±0.128 
1.268 ±0.152 
1.018 ±0.097 
0.947 ±0.069 
0.844 ± 0.088 
0.814 ± 0.048 



To correctly determine the efficiency and acceptance using the MC simulation, 
the kinematic distributions in the MC simulation should match the data. As an 
example Figure [4] shows a comparison between data and MC of the muon for the 
D°(K~ yU + z/)7T + analysis in each trigger period. 

Since there is no clear evidence of D° — > fi + n~ an upper limit on N(fifi) is deter- 
mined assuming that the number of events found in the signal region is the sum of 
signal and background events (both obeying Poisson statistics). 

The hadronic decays D° — > K~ir + and D° — > 7r + 7r~ can result in an extra con- 
tribution to the low mass sideband and signal region, respectively. This potential 
contamination was measured in data using reconstructed D° — > K~tt + decays, the 
measured CMS misidentification rate ~ 0.01% [9], and, for the D° — > 7r + 7r~ decay, the 
known relative branching ratio [1]. It is found that the hadronic D° decays produce 
a negligible contribution to the signal and sideband regions. 

From Figure [2] it is clear that the background in this limited mass region can be 
assumed to be flat, and an estimate can be made from the sidebands. Three different 
regions are defined corresponding to the signal region, M(D°)pdg ± 40 MeV, where 
40 MeV corresponds to ~ 2.1 times the mass resolution measured in MC, and two 
sidebands (40 MeV wide). The ranges of these three regions are: 1.825 < M MM < 1.905 
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Figure 4: Comparison between data and MC of the muon transverse momentum for 
the D°(K~ yU + z/)7r + analysis in each trigger period. 
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for the signal region, 1.765 < M w < 1.805 for the left sideband and 1.925 < M w < 
1.965 for the right sideband. We find 23 background events (A^g) and 23 events 

observed in the signal region (A^g). In more detail the D° — > K~ n + v yields and the 
number of n + fi~ events in the signal and sideband regions are shown for each period 
in Table H 



Table 2: Yields for D° — > K jj, + iy, N ^ s and A^g for D° — > fi + /i , and systematic 
uncertainty in the seven periods. 



Trigger 




^obs' ^bkg 


Systematic uncertainty 


HLT_Mu3 


2412 ± 145 


0, 


19.1% 


HLT_Mu5 


2447 ± 357 


1, 


18.0% 


HLT_Mu7 


11799 ±215 


6, 4 


19.0% 


HLT_Mu9 


9982 ± 176 


6, 6 


17.4% 


HLTJvlull 


10079 ± 185 


3, 5 


16.4% 


HLT_Mul5 (2010) 


5302 ± 118 


1, 3 


17.9% 


HLT_Mul5 (2011) 


16458 ± 204 


6, 5 


15.6% 



The limit on the D° — > jj, + [i~ branching fraction depends on statistical and sys- 
tematic uncertainties. Several sources of systematic uncertainties have been con- 
sidered. There is no systematic uncertainty on the production cross section, as we 
use the ratio: ^ t 2~*Sn& The statistical errors on the determination of the 

D*+— >D U (K~ ij 1 +i/)it+ 

acceptance and of the reconstruction efficiency ratio e Tec (K fxu) / e Tec (fj,fj,) (see Table 
[T| are taken as a systematic uncertainties. The uncertainties on the reconstruction 
efficiency of kaons and muons have been determined using data and found to be 
3.9% for the hadron tracking efficiency and 1% for the muon tracking efficiency [10]. 
The biggest systematic uncertainty comes from the determination of the trigger ef- 
ficiency, e^rig. ^he reason is that the Pt spectrum of muons coming from charm 
meson decays is not matched by the CMS triggers. The consequence is an extremely 
low e^ r j g , especially for the normalization decay mode, which can be a significant 
source of systematic uncertainty. To estimate this contribution we compute the ratio 
-^trig = (Yieldrj a t a /Luminosity)/et t for the different triggers. If the Monte Carlo 
correctly simulates the trigger, the ratio -R^rig should not depend on the different 
HLT triggers. The weighted average of the seven periods is calculated as well as the 
PDG scale factor S [4J, S = ^/(DOF - 1) = 1.131. This S-factor tells how much 
any single error should be increased to have xV(DOF — 1) = 1. We estimate an 
uncertainty of 13.1% for this source. Variations of the fitting functions used to obtain 
the yield of the normalization mode gives systematic uncertainties of 1% to 9% for 
the seven run periods. Another source of systematic uncertainties is the contami- 
nation from other decay modes to the Yield(D* + — > D° (K~ [i + v)-k + ) . We consider 



D° — > K*~ (K~7r )fi + v which is the largest contamination for the semileptonic decay 
D° — > K~fj, + v. Using a MC simulation we estimate a contribution of 1.8% from 
D° — > K*~ \K~ it ) /i + v . Finally, the uncertainty on the PDG value of the branching 
fraction B(D° — > K~ fi + v) is included. Adding these contributions in quadrature we 
obtain the systematic uncertainty for each period shown in Table |2j which has been 
included in the determination of the upper limit. 



3 Conclusions 



The 90% confidence level upper limit is computed using the CL S approach [HI [12], 
combining the results of the 7 periods. The values used are those reported in Table 
|2jas well as e^g(K fiu) / e^g(fifi) and erec(-^A tz/ )/ e rec(A i A i ) as shown in Table 
The final result is: 



B(D° //V~) < 5.4 x 1(T 7 (90% CL). (2) 

In summary, the FCNC decay D° — > n + fi~ has been searched for using the CMS 
detector. No evidence has been found in ~ 90pb _1 of data. We show the upper limit 
at 90% confidence level in Table |3j together with the present published best limits. 
Although this upper limit is not the best limit for this FCNC decay, it is the first 
time a semileptonic decay has been used as the normalization. 



Table 3: Upper limit at 90% CL for D° -> /i+/i~. 



Experiment 


Upper limit at 90% CL 


BABAR [13] 
CDF [H] 
BELLE [15] 
this measurement 


< 1.3 x 10~ 6 

< 2.1 x 10~ 7 

< 1.4 x 10~ 7 

< 5.4 x 10~ 7 
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